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Partially  shaded  conditions  (PSCs)  often  occur  in  large  photovoltaic  generation  systems  (PGSs).  PSCs 
cause  losses  in  system  output  power,  hot  spot  effects,  and  system  safety  and  reliability  problems.  When 
PSC  occur,  the  PGS  power- voltage  characteristic  curve  exhibits  multiple  peak  values;  that  is,  the  curve 
comprises  a  global  maximum  power  point  and  multiple  local  maximum  power  points.  Current  literature 
includes  various  studies  of  global  maximum  power  point  tracking  (GMPPT)  algorithms  and  hardware 
architectures  suitable  for  PSC;  because  the  substantial  quantity  of  PSC  literature,  this  subject  must  be 
comprehensively  reviewed.  To  focus  on  GMPPT  techniques  used  in  PSC,  traditional  maximum  power 
point  tracking  techniques  and  circuit  architectures  that  cannot  distinguish  GMPP  and  LMPP  were  not 
discussed. 
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1.  Introduction 

According  to  a  technical  report  by  the  European  Photovoltaic 
Industry  Association  (EPIA),  more  than  31  GW  of  photovoltaic  (PV) 
capacity  was  installed  globally  in  2012;  aggressive  projections 
indicate  that  this  number  could  increase  to  84  GW  in  2017  1  . 
Weather  and  environmental  factors  affect  the  output  character¬ 
istics  of  PV  cells,  that  is,  cell  output  voltage  and  current  vary 
according  to  changes  in  irradiance  and  temperature.  Therefore,  a 
unique  maximum  power  point  (MPP)  exhibits  under  specific 
irradiance  and  temperature.  The  maximum  power  point  tracking 
(MPPT)  algorithm  is  crucial  in  attaining  the  maximal  PV  power, 
facilitating  optimal  PV  cell  performance.  Researchers  have  pro¬ 
posed  numerous  MPPT  algorithms  that  demonstrate  excellent 
tracking  efficiency  in  uniform  insolation  conditions.  However,  PV 
generation  systems  (PGSs)  in  urban  environments  installed  on 
rooftops  or  eaves  are  easily  shaded  by  neighboring  buildings, 
clouds,  or  dust  coverage,  causing  PGSs  to  generate  power  output 
that  differs  from  that  in  ideal  insolation  conditions.  Shading  is 
especially  problematic  regarding  PGSs  that  exhibit  wide-range 
layouts  primarily  because  complete  shading  of  a  large  area  is 
unlikely,  causing  partially  shaded  conditions  (PSC).  When  PSC 
occur,  certain  PV  cells  are  exposed  to  normal,  non-shaded  condi¬ 
tions,  and  output  normal  power  levels;  however,  other  cells  yield 
decreased  output  levels  because  of  shading,  generating  a  mis¬ 
match  in  the  overall  PGS.  Because  of  the  characteristics  of  PV  cells, 
shaded  conditions  significantly  reduce  cell  output  current  but  not 
output  voltage.  This  mismatch  in  power  generation  influences 
series-connected  PGSs;  however,  traditional  MPPT  algorithms  are 
designed  to  address  only  uniform  insolation  and  cannot  ade¬ 
quately  account  for  PSC.  Consequently,  developing  new  MPPT 
techniques  suitable  for  application  in  PSC  is  critical.  Numerous 
studies  regarding  this  subject  have  been  published;  thus,  compos¬ 
ing  an  accurate  and  comprehensive  reference  document  that 
details  information  related  to  various  MPPT  techniques  applicable 
to  PSC  is  essential. 

Various  reviews  [2-10]  of  MPPT  methods  for  PGSs  exist  in  the 
current  literature;  of  these  studies,  most  of  the  reviews  [2-7] 
researched  only  traditional  MPPT  methods,  and  Ishaque  and  Salam 
[8]  and  Salam  et  al.  [9]  focused  on  both  traditional  methods  and 
those  suitable  for  PSC.  When  PSC  occur,  the  power-voltage  (P-V) 
characteristic  curve  of  the  PV  modules  becomes  complex,  exhibit¬ 
ing  multiple  peak  values.  Therefore,  traditional  MPPT  methods  are 
not  suitable  for  use  in  PSC.  To  provide  PGS  researchers  and 
engineers  with  useful  reference  information,  this  study  explored 
global  maximum  power  point  tracking  (GMPPT)  methods  suitable 
for  application  in  PSC.  First,  the  basic  operating  principles  of 
various  GMPPT  methods,  including  their  advantages  and  limita¬ 
tions,  were  described.  Regarding  firmware-related  techniques, 
various  indices  were  evaluated,  including  whether  the  design 
parameters  were  relevant  to  the  system,  GMPPT  time,  GMPPT 
accuracy,  algorithm  complexity  and  required  sensors.  In  terms  of 
hardware-related  techniques,  the  advantages  and  disadvantages  of 
various  architectures  were  discussed.  Finally,  an  in-depth  simula¬ 
tion  on  some  highly  cited  GMPPT  methods  was  made  and  conclu¬ 
sion  was  offered  accordingly.  It  was  difficult  to  propose  a  fair 
benchmark  for  various  techniques  mentioned  in  the  literature 
primarily  because  researchers  used  distinct  PGSs  (including  power 
rating  and  power  converter  topology)  and  shading  patterns  to  test 
their  proposed  algorithms.  Nonetheless,  the  authors  of  the  current 
study  endeavored  to  employ  several  fair  evaluation  indices  and 
verify  the  correctness  of  this  performance  evaluation  using  the 
simulated  results  in  Section  4.  The  information  gathered  herein 
should  provide  researchers  with  comprehensive  reference  infor¬ 
mation  and  recommendations  for  future  studies  addressing  PSC. 
Although  Bidram  et  al.  [10]  also  compared  GMPPT  methods,  this 


study  additionally  classified  various  GMPPT  methods  into  cate¬ 
gories  and  explored  their  usage  limitations  and  design  considera¬ 
tions.  Furthermore,  shading  detection  methods  proposed  in  the 
literature  were  compared,  and  a  simulation  platform  was  provided 
for  verifying  the  performance  of  some  highly  cited  GMPPT  meth¬ 
ods.  To  the  best  of  the  authors’  knowledge,  no  existing  literature 
has  provided  this  type  of  categorizing  and  testing. 

2.  Exploration  of  the  effects  of  partially  shaded  conditions 

When  PV  cells  are  connected  in  series,  a  bypass  diode  is 
typically  added  in  parallel  for  protection,  preventing  damage  due 
to  overheating  when  one  or  more  PV  cells  are  shaded  or  shorted. 
However,  adding  a  bypass  diode  changes  the  current-voltage  (/-V) 
characteristic  curve  when  the  PGS  is  shaded  by  environmental 
conditions  such  as  clouds,  tree  shade,  and  building  shadows. 
Furthermore,  the  P-V  characteristic  curve  also  changes  accordingly 
and  exhibits  multiple  peak  values.  Current  studies  on  the  effects  of 
PSC  can  be  divided  into  three  categories:  (1)  deriving  analytical 
mathematical  formulas  to  describe  the  I-V  and  P-V  characteristic 
curve  [11,12];  (2)  establishing  a  simulation  or  experimental  plat¬ 
form  that  can  evaluate  the  effects  of  PSC  [13-20];  and  (3)  evaluat¬ 
ing  the  power  losses  caused  by  PSC  [21-25  .  This  section  explores 
the  effects  of  PSC  in  PGSs. 

2.1.  Basic  characteristics  of  PV  cells 

Based  on  their  electrical  characteristics,  PV  cells  are  equivalent 
to  a  one-diode  model  as  shown  in  Fig.  1.  This  model  includes  a 
current  source  7S,  a  diode  current  7d,  an  equivalent  parallel 
resistance  RP,  and  an  equivalent  series  resistance  Rs.  The  output 
current  of  the  PV  module  IPV  is  the  difference  between  the 
photogenerated  current  7S  and  diode  current  7d,  as  shown  in  (1): 

w.,s-,„(exp(S „> 

where  7So  is  the  saturation  current,  VPV  is  the  output  voltage  of  PV 
module,  q  is  the  charge  of  a  single  electron,  a  is  the  ideal  factor  of 
the  diode,  k  is  the  Boltzmann  constant,  T  is  the  temperature  in 
Kelvin. 

2.2.  Effect  of  partial  shade  conditions  on  photovoltaic  generation 
systems 

PGSs  are  formed  by  several  PV  modules,  which  comprise 
several  PV  cells  connected  in  series  or  parallel.  Parallel  connections 
can  increase  the  output  current  of  the  module,  and  serial  connec¬ 
tions  can  improve  the  output  voltage.  When  PSC  occur,  the  shaded 
PV  cell  becomes  a  load  instead  of  a  generator,  producing  a  hot  spot 
that  can  subsequently  destroy  PV  cells.  Therefore,  commercially 
available  PV  modules  often  comprise  a  parallel-connected  bypass 
diode  to  prevent  this  problem.  When  PSC  occur,  currents  flow 
through  the  bypass  diode  instead  of  the  shaded  PV  module, 
generating  multiple  steps  in  the  7-V  characteristic  curve  and 


Fig.  1.  Equivalent  circuit  of  photovoltaic  cells. 
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Pattern  1 


Pattern  2 


Pattern  3 


Fig.  2.  A  5slp  PGS  with  various  shading  patterns. 


a  b 


Fig.  3.  The  characteristic  curves  of  the  5slp  PGS  in  Fig.  2:  (a)  P-V  characteristic  curve;  (b)  I-V  characteristic  curve. 


multiple  peak  values  in  the  P-V  characteristic  curve.  For  example, 
Fig.  2  shows  a  PGS  comprising  five  PV  modules  connected  in  series 
(denoted  as  5slp  system).  Its  P-V  and  I-V  characteristic  curves  for 
different  shading  patterns  are  illustrated  in  Fig.  3(a)  and  (b), 
respectively.  Fig.  3  shows  that  GMPP  under  PSC  may  occur  on  a 
larger  voltage  range  compared  with  that  in  PGSs  under  uniform 
insolation. 


3.  Global  maximum  power  point  tracking  techniques  for 
partial  shading  conditions 

In  this  section,  GMPPT  techniques  suitable  for  PSC  are  explored. 
The  solutions  for  PSC  proposed  in  current  literature  can  be  divided 
into  firmware-based  [26-69],  and  hardware  architecture-based 
solutions  [70-113].  Fig.  4  displays  common  PGS  architectures, 
illustrating  that  when  PGSs  use  a  centralized  architecture,  only 
one  set  of  power  converter  is  required;  therefore,  firmware  must 
be  used  to  address  the  multiple  peak  values  that  the  PSC  cause  in 
P-V  characteristic  curves.  On  the  other  hand,  when  PGSs  use 
distributed  architectures,  multiple  sets  of  power  converter  exist; 
therefore,  each  power  converter  can  track  its  own  MPPT.  Because 


hardware-based  solutions  relate  to  the  power  converter  topo¬ 
logy  and  system  design,  this  study  mainly  focuses  on  exploring 
firmware-based  GMPPT  techniques.  Fig.  5  shows  the  P-V  charac¬ 
teristic  curve  of  a  typical  5slp  system  under  a  specific  shading 
pattern;  this  system  was  used  as  the  basis  to  clarify  operating 
concepts  and  design  guidelines  regarding  various  GMPPT  algo¬ 
rithms.  (Fig.  6). 

3.2.  Firmware-based  GMPPT  techniques 
3.1.1 .  Two-stage  method 

Fig.  5  shows  that  when  PSC  occur,  the  P-V  characteristic  curve 
of  PGS  exhibits  multiple  peak  values,  and  the  number  of  peak 
values  correlates  with  the  shading  pattern.  By  dividing  this  P-V 
characteristic  curve  into  multiple  single-peak  P-V  characteristic 
curves,  traditional  MPPT  methods  (e.g.,  perturb  &  observe  (P&O)  or 
incremental  conductance  (IncCon)  methods)  can  be  used  to 
determine  the  maximal  value  of  each  single-peak  curve.  Therefore, 
the  two-stage  technique  involves  using  methods  proposed  in  the 
first  stage  to  determine  the  approximate  GMPP  location  (i.e.,  an 
interval),  subsequently  using  traditional  MPPT  methods  to  deter¬ 
mine  the  precise  GMPP  location.  Thus,  the  key  to  the  success  of  the 
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Fig.  4.  Common  PGS  architectures. 


Fig.  5.  P-V  characteristic  curve  under  a  specific  shading  pattern. 


two-stage  method  is  the  accuracy  of  the  GMPP  interval  obtained  in 
the  first  stage.  The  following  sections  clarify  the  two-stage  method 
used  in  the  literature. 

3J.2.2.  System  characteristic  curve  method.  Works  in  [26-28] 
employed  a  similar  technique,  using  a  preset  linear  function  to 
move  the  operating  point  (OP)  near  the  GMPP  under  PSC,  as 
shown  in  Fig.  7.  This  linear  function  is  related  to  various  system 
parameters,  including  open  circuit  voltage  and  short  circuit 
current.  The  linear  function  used  in  [26]  can  be  expressed  as 
follows: 

n  _  Vpm 

Kpm  —  j  \2) 

ipm 
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Fig.  6.  Concept  of  dividing  multiple-peak  P-V  characteristic  curve  into  several 
single-peak  curve. 

The  linear  function  used  in  [27]  can  be  expressed  as  follows: 
PMPp(t)  =  a[T(t)]lMPp[E(t)]  (3) 

And  the  linear  function  used  in  28]  can  be  expressed  as  follows: 

V*pv  =  fef)  x  Ipv[n]  *  (p)  x  Ipv[n]  (4) 

The  advantage  of  this  method  is  its  tracking  speed.  Because  the 
first  stage  only  requires  one  step  to  move  the  OP  to  the  vicinity  of 
the  GMPP,  the  tracking  time  of  this  method  is  similar  to  that  of 
traditional  MPPT  methods.  However,  disadvantages  exist.  First, 
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Fig.  7.  System  characteristic  curve  method. 
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Fig.  8.  Incorrect  tracking  of  GMPP  using  the  system  characteristic  curve  method. 


Fig.  10.  Incorrect  tracking  of  GMPP  using  the  method  proposed  by  Patel. 


obtaining  the  open  circuit  voltage  and  short  circuit  current  of  the 
system  requires  open  or  short  circuits.  Implementing  this  in  an 
actual  system  can  cause  power  losses  or  safety  concerns.  Second, 
although  the  linear  function  used  in  this  method  can  move  the  OP 
near  the  GMPP  in  certain  shading  patterns,  an  erroneous  location 
of  the  new  OP  can  prevent  tracking  the  GMPP;  thus,  this  method 
cannot  guarantee  that  the  GMPP  can  be  tracked.  This  concept  is 


shown  in  Fig.  8.  Furthermore,  when  shading  patterns  become 
complex,  the  probability  that  this  method  can  successfully  track 
GMPP  decreases.  Ref.  [29  is  one  conference  paper  that  used  a 
concept  similar  to  this  method. 


3. 2. 2.2.  Two-stage  searching  method.  In  a  two-stage  searching 
method,  the  first  stage  is  used  to  search  for  the  GMPP  location 
interval;  traditional  MPPT  methods  are  then  used  in  the  second 
stage  to  find  the  precise  GMPP  location.  By  using  this  method, 
Patel  and  Agarwal  [30]  employed  various  shading  pattern 
observations  as  clues  to  determine  the  GMPP  location.  As  shown 
in  Fig.  9,  the  basic  search  rules  for  [30]  are  as  follows: 

1 )  0.85  Voc,aii  (Voc,aii  represents  the  total  open  circuit  voltage  of  the 
system)  is  used  as  the  P&O  search  starting  point  to  find  and 
record  a  peak  value. 

2)  The  peak  value  found  in  Step  1  is  used  as  the  basis  for  moving 
the  OP  one  large  step  (the  reference  recommended  0.6-0.7  Voc, 
one  (Voc, one  represents  the  open  circuit  voltage  of  a  single 
module))  to  the  left,  then  conduct  P&O  to  search  and  record 
the  next  peak  value. 

3)  If  the  obtained  peak  value  is  higher  than  the  previous  peak 
value,  Step  2  is  repeated;  if  the  determined  peak  value  is 
smaller  than  the  previous  peak  value,  then  the  previous  peak 
value  is  the  GMPP. 

Ref.  [30]  is  the  paper  with  the  highest  number  of  citations 
among  similar  literature;  its  method  is  easy  to  implement,  and  can 
be  integrated  into  traditional  PGS  firmware;  however,  disadvan¬ 
tages  exist.  First,  the  authors  of  [30^  found  that  when  P-V  curve  is 
traversed  from  either  side,  the  magnitude  of  the  peaks  increases; 
after  reaching  the  GMPP,  the  magnitude  of  the  subsequent  peaks 
continuously  decreases.  Therefore,  Step  3  can  be  used  to  deter¬ 
mine  the  GMPP.  However,  the  current  study  explored  all  possible 
shading  patterns  and  found  that  there  exists  some  exception  about 
this  claim  (for  detailed  descriptions,  please  refer  to  Section  4). 
Consequently,  this  method  cannot  ensure  successful  GMPPT. 
Fig.  10  shows  one  example  of  this  exception.  In  addition,  the 
parameter  used  in  this  method  (i.e.  value  of  the  large  step)  affects 
the  probability  of  tracking  the  GMPP.  Another  disadvantage  of  this 
method  is  the  tracking  speed.  Because  each  LMPP  must  be 
determined  using  P&O  method,  this  method  requires  more  track¬ 
ing  steps  compared  with  other  methods. 

Works  in  [31-33]  used  similar  methods,  employing  a  large 
interval  to  search  the  entire  P-V  characteristic  curve  (called 
global-biased  search)  to  determine  the  largest  peak  value.  Next, 
a  refined  search  is  conducted  around  the  largest  peak  value  found 
in  previous  stage  to  locate  the  GMPP  (called  local-biased  search). 
Of  these  three  approaches,  Renaudineau  et  al.  [31]  used  a, 
Koutroulis  et  al.  [32]  used  fixed  power  intervals,  and  Kouchaki 
et  al.  [33]  used  0.8  Voc,one  as  the  step  for  global-biased  search.  The 
concept  of  this  kind  of  method  (Ref.  [32]  by  Koutroulis)  is  shown  in 
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Fig.  11.  Method  proposed  by  Koutroulis. 


Fig.  12.  Incorrect  tracking  of  GMPP  using  the  method  proposed  by  Koutroulis. 


Fig.  11.  The  advantage  of  this  kind  of  method  is  its  simplicity;  it  can 
easily  be  integrated  with  traditional  PGS  firmware  (requiring  only 
an  additional  global-biased  search  stage).  The  disadvantage  of  this 
kind  of  method  is  that  the  selected  global-biased  search  step 
affects  the  probability  of  finding  the  GMPP.  When  the  selected  step 
is  large,  the  system  requires  comparatively  less  time  to  locate  the 
GMPP,  but  may  miss  the  GMPP,  as  shown  in  Fig.  12.  Conversely, 
when  a  small  step  is  selected,  the  system  has  an  increased 
probability  of  finding  the  GMPP,  but  requires  relatively  more  time. 
To  use  this  method,  a  tradeoff  between  the  search  step  size  and 
the  time  required  must  be  addressed  (for  details,  please  refer  to 
Section  4).  In  addition,  the  method  in  [32]  used  power  as  the  step 
variable  for  the  search.  Although  this  thinking  is  intuitive  and 
yields  a  comparatively  fast  tracking  speed,  this  method  requires 
using  a  specifically  designed  power  converter  that  can  accept 
power  command  for  constant  power  control.  Kashif  et  al.  34]  and 
Kouchaki  et  al.  [35]  are  conference  papers  used  similar  two-stage 
searching  concepts. 


3.2.2.  Segmental  search  method 

The  basic  principle  of  the  segmental  search  method  is  similar  to 
that  of  the  two-stage  searching  method  described  in  3.I.I.2.  The 
primary  difference  is  that  two-stage  searching  method  uses 
traditional  MPPT  methods  (such  as  P&O)  as  the  basis  for  the 
second  stage  search.  Segmental  search,  however,  is  based  on 
mathematical  theories,  finding  the  GMPP  by  gradually  reducing 
the  search  range. 


3. 2.2. 2.  DIRECT  method.  The  theoretical  basis  for  the  method  used 
in  [36  is  that  the  P-V  characteristic  curve  of  PV  modules  conforms 
to  Lipschitz  characteristics;  therefore,  the  dividing  rectangle 
(DIRECT)  technique  can  be  used  to  find  the  peak  values.  DIRECT 
refers  to  dividing  the  search  area  into  three  areas  of  equal  intervals 
by  using  mathematical  methods,  subsequently  using  mathematical 
equations  to  find  potentially  optimal  intervals  (POIs)  that  exhibit 


potential  peak  values.  The  rules  for  POI  search  are  as  follows: 

/(Xj)  +  K^LZ^>/(Xi)  +  K^=^  Vl 


(aj  -  bj) 


—  f  max  ~P£\f 


max 


(5) 

(6) 


The  search  area  is  divided  into  three  equal  portions  to  locate 
the  POIs  before  further  subdivision,  as  shown  in  Fig.  13. 

Ref.  [36  has  been  frequently  cited  and  compared.  The  advan¬ 
tage  of  this  method  is  that  it  is  based  on  a  solid  mathematical 
foundation  and  facilitates  rapid  tracking;  its  disadvantages  are 
similar  to  those  of  the  two-stage  searching  method.  In  certain 
shading  patterns,  this  method  cannot  track  the  GMPP,  although 
this  is  unlikely.  Another  potential  disadvantage  is  that  although 
the  mathematical  equations  used  in  this  method  are  relatively 
simple  and  the  segmentation  concept  can  be  easily  implemented, 
the  author  added  numerous  exceptions  during  actual  implemen¬ 
tation  to  avoid  missing  the  GMPP.  For  example,  the  author 
proposed  conducting  nine  interval  searches  within  the  search 
range;  however,  if  the  POIs  that  contained  the  two  largest  values 
were  not  adjacent  to  each  other,  then  a  wide  area  search  with  27 
intervals  had  to  be  re-conducted.  Conversely,  if  the  obtained  POIs 
were  adjacent  to  each  other,  then  segmentation  could  be  con¬ 
ducted  on  the  POI  with  the  maximal  value.  This  type  of  exception 
adds  complexity  and  difficulty  when  firmware  engineers  imple¬ 
ment  the  algorithm  based  on  the  DIRECT  method.  In  addition, 
this  method  cannot  be  directly  integrated  into  traditional  PGS 
firmware. 


3.122.  Fibonacci  methods.  Works  in  [37]  and  [38]  used  Fibonacci 
sequences  as  a  basis  for  GMPP  searches.  Compared  with  the 
DIRECT  method  in  3.1.2.1,  the  primary  difference  of  this  method 
is  using  the  Fibonacci  sequence  as  the  mathematical  basis  for 
segmentation,  as  shown  in  Fig.  14.  The  equation  is  as  follows: 


*2  =*3  + 


’F(n-l)" 
.  F(n)  . 


(X4-X3) 


Fig.  13.  The  concept  of  DIRECT  method. 
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Fig.  14.  The  concept  of  the  Fibonacci  method. 
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Fig.  15.  ANN  with  three  layers  (one  hidden  layer). 


X\  =  X4  + 


T(n- 1)' 
.  F(n)  _ 


(X4-X3) 


expressed  as  follows: 

5 

n]  =  2  WjjUf  +  bj 

i  =  1 


where  F(n- 1)  and  F(n)  represent  the  (n-  7  )th  and  nth  number 
in  the  Fibonacci  sequence.  The  advantages  and  disadvantages  of 
this  method  are  the  same  as  those  in  the  DIRECT  method. 

Ref.  [39]  and  Ref.  [40]  are  conference  papers  related  to  the 
segmental  search  method.  Fan  et  al.  [39]  used  a  random  number 
method  to  randomly  sample  power  values  under  six  distinct 
voltages,  subsequently  using  the  point  with  the  highest  power 
value  (of  the  six  selected  points)  to  determine  the  new  search 
range.  Repeat  random  sampling  was  conducted  on  the  new  range 
until  convergence  was  achieved.  Escobar  et  al.  [40]  used  the 
geometric  characteristics  of  the  I-V  characteristic  curve  to  conduct 
a  sweep  search  with  variable  steps. 


3.1.3.  Soft  computing  method 

Soft  computing  (SC)  refers  to  using  computers  to  simulate  the 
biochemical  processes  of  natural  intelligence  systems  (e.g.,  human 
perception,  brain  structure,  and  evolutionary  process)  to  solve 
non-deterministic  polynomial-time  hard  (NP-hard)  or  NP- 
complete  problems.  Because  finding  the  GMPP  in  multi-peak  P-V 
characteristic  curves  caused  by  PSC  can  be  viewed  as  an  optimiza¬ 
tion  problem,  scholars  have  used  SC  techniques  to  address  such 
problems.  SC  methods  include  artificial  neural  networks  (ANNs), 
fuzzy  logic  control  (FLC),  evolutionary  computation  (including 
genetic  algorithms  (GAs),  differential  evolution  (DE),  particle 
swarm  optimization  (PSO),  ant  colony  systems  (ACSs)),  and  chaotic 
search  (CS). 


3. 2.3. 2.  Artificial  neural  network.  An  ANN  is  a  computational  model 
that  uses  interconnected  artificial  neurons  to  imitate  the  ability  of 
biological  neural  networks.  An  ANN  with  three  layers  is  shown  in 
Fig.  15.  From  Fig.  15,  users  can  input  information  (such  as 
irradiance  level  and  the  temperature  of  PV  modules)  into  the 
input  layer  and  calculate  the  output  (such  as  the  GMPP  location) 
according  to  the  weighting  values  of  ANN.  For  example,  the  input 
and  output  relation  of  an  ANN  as  shown  in  Fig.  15  can  be  obtained 
using  (9)  to  (12).  The  neurons  in  the  input  layer  are  responsible 
only  for  transmitting  the  input  information  to  the  hidden  layer. 
The  input  of  the  number  j  unit  in  the  hidden  layer  can  be 


where  wjt  is  the  weighting  of  the  ith  input  unit  connected  to  the 
jth  hidden  layer  unit,  and  bj  represents  the  bias  value  of  the  jth 
unit  in  the  hidden  layer.  The  output  of  the  jth  unit  in  the  hidden 
layer  can  be  expressed  as  follows: 

Vj=f,<nJ)  (10) 

where  fj  is  the  activation  function  used  in  this  layer.  Commonly 
used  activation  functions  include  liner  function,  hyperbolic  tan¬ 
gent  sigmoid,  and  Log-sigmoid.  Notably,  if  a  low-cost  microcon¬ 
troller  with  low  calculation  capability  is  used  to  implement  the 
ANN,  the  linear  function  is  preferred  because  it  is  relatively  easy 
to  apply. 

The  input  of  the  kth  unit  in  the  output  layer  can  be  expressed  as 
follows: 

9  10  2 

K=  E  wkjUj  +  bk  (11) 

2  =  1 

where  wkj  is  the  weighting  value  of  the  jth  hidden  layer  unit 
connected  to  the  kth  output  unit,  and  bfc  represents  the  bias  value 
of  the  kth  output  unit.  The  output  of  the  kth  output  unit  can  be 
expressed  as  follows: 

Vk=h(n2k)  (12) 

Similarly,  f2  is  the  activation  function  used  for  the  output  layer. 

By  using  large  amounts  of  training  data,  the  ANN  continually 
adjusts  the  weighting  and  bias  values,  allowing  the  network- 
calculated  output  to  approximate  the  target  output.  One  of  the 
most  common  ANN  learning  methods  is  the  back-propagation 
method.  The  advantage  of  ANNs  is  their  parallel  computing 
capability.  Other  SC  methods  may  require  multiple  iterations  to 
obtain  the  optimal  solution,  whereas  ANNs  can  use  simple  multi¬ 
plication  and  addition  to  rapidly  calculate  output.  Therefore,  ANNs 
enable  rapid  calculation.  However,  the  accuracy  of  an  ANN  is 
determined  based  on  its  training  data.  If  the  training  data  are 
insufficient,  or  the  data  do  not  cover  the  entire  problem  space, 
then  the  accuracy  of  the  ANN  will  be  reduced  accordingly. 
Currently,  the  ANN-based  MPPT  method  can  only  be  used  in 
uniform  insolation  conditions  [41,42]  primarily  because  at  a  single 
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irradiance  level,  the  location  of  the  MPP  is  only  related  to  the 
irradiance  and  the  temperature.  Thus,  the  training  data  are  simple. 
When  PSC  occur,  the  irradiance  of  the  module,  module  tempera¬ 
ture  and  shading  pattern  all  affect  the  MPP  location.  Consequently, 
the  training  data  required  by  the  ANN  substantially  increase,  and 
these  data  are  not  easily  collected.  This  is  the  primary  limitation  of 
applying  ANNs  to  GMPP  searches.  Therefore,  only  works  in  [43] 
and  [44]  in  current  literature  proposed  using  ANNs  as  the  primary 
method  for  addressing  PSC.  The  input  variables  used  in  [43]  and 
[44]  were  the  averaged  irradiance  of  selected  modules,  and  the 
design  was  affected  by  the  arrangement  method  of  PV  modules. 
Therefore,  if  the  architecture  of  the  PV  array  changes  (e.g.,  adding 
new  PV  panels),  the  ANN  must  be  re-trained.  In  addition,  irra¬ 
diance  and  temperature  sensors  are  more  expensive  compared 
with  the  voltage  and  current  sensors  used  in  other  GMPP  methods. 
Ref.  [45]  is  one  conference  paper  that  used  ANN  concepts. 

3 A3. 2.  Fuzzy  logic  control  Traditional  control  system  design 
requires  understanding  the  system  being  controlled;  that  is,  using 
precise  mathematical  models  to  describe  the  system.  However,  when 
the  controlled  system  becomes  overly  complex,  it  is  often  difficult  to 
use  system  identification  method  to  establish  a  system  model.  On  the 
other  hand,  FLC  uses  fuzzy  set  theory  to  convert  the  linguistic  values 
directly  into  automatic  control  strategies  using  expert  knowledge  or 
operator  experience.  The  advantages  of  FLC  include  the  following:  1 ) 
no  precise  mathematical  model  is  required;  and  2)  FLC  can  integrate 
the  knowledge  of  human  experts  into  the  controller  design  process. 
When  conducting  fuzzy  control,  the  behavior  of  the  system  being 
controlled  is  described  by  a  set  of  fuzzy  rules.  These  fuzzy  rules  use 
semantic  fuzzy  information  rather  than  mathematical  equations. 
Therefore,  FLC  can  convert  the  knowledge  of  human  experts  into 
fuzzy  control  rules,  reducing  the  complexity  of  control  system  design. 
Fig.  16  shows  a  typical  block  diagram  of  FLC.  From  Fig.  16,  a 
FLC  comprises  three  primary  parts:  fuzzification  stage,  a  fuzzy 
inference  engine,  and  defuzzification  stage.  In  the  fuzzification 
stage,  membership  function  is  used  to  convert  the  input  data  into 
linguistic  values.  The  fuzzy  inference  engine  must  be  used  in 
conjunction  with  the  fuzzy  rule  base,  and  can  be  used  to  calculate 
the  required  control  quantity.  In  the  defuzzification  stage, 
mathematical  equations  are  used  to  convert  the  output  variable 
from  a  linguistic  value  to  a  crisp  value.  The  center  of  gravity  and 
center  of  area  methods  are  often  used  for  defuzzification.  Regarding 
FLC,  establishing  a  rule  base  is  critical  in  the  design.  The  if-then  rules 
formed  using  linguistic  variables  can  be  employed  to  establish  the 
FLC  system  input  and  output  relations,  and  must  be  designed  to 
agree  with  the  knowledge  of  experts  concerning  a  specific  problem. 
FLC  is  particularly  suitable  for  use  in  non-linear,  time-varying 
systems  and  systems  with  incomplete  models.  Thus,  FLC  can  be 
used  to  deal  with  PSC  problems.  Karatepe  et  al.  [46]  used  distributed 
architecture  and  utilized  FLC  to  replace  traditional  MPPT  methods. 


Because  each  power  converter  contains  a  designated  MPPT  controller, 
the  system  guarantees  that  the  GMPP  can  be  tracked.  The  advantage 
of  the  method  used  in  [46]  is  its  rapid  tracking  speed  and  high 
tracking  accuracy;  its  disadvantage  is  that  it  exhibits  a  higher 
hardware  cost  compared  with  that  of  centralized  architectures. 
Alajmi  et  al.  [47]  used  a  method  similar  to  the  two-stage  searching 
method  in  3.1.1.2  to  conduct  GMPPT.  This  method  first  swept  the  P-V 
characteristic  curve  and  recorded  various  LMPPs,  then  replaced  P&O 
by  using  FLC  to  conduct  GMPPT.  This  method  yields  a  rapid  tracking 
speed,  but  demonstrates  disadvantages  similar  to  those  of  other  two- 
stage  searching  methods.  Ref.  [48  is  one  conference  paper  that  used 
FLC  concepts. 

3 A. 3. 3.  Genetic  algorithm.  A  GA  is  an  evolutionary  computation 
(EC)  method  that  mimics  the  “survival  of  the  fittest”  evolutionary 
law  in  the  natural  world;  its  primary  operators  are  reproduction, 
crossover,  and  mutation.  The  process  of  using  GA  to  obtain  the 
optimal  solution  for  a  problem  is  as  follows:  1)  the  search 
parameters  are  coded  as  binary  strings  in  a  chromosome;  2)  this 
process  is  randomly  repeated  to  produce  N  initial  species  (strings); 
3)  the  fitness  function  is  designed  based  on  the  solution 
conditions.  Species  that  exhibit  high  fitness  values  are  selected 
for  the  mating  pool  (i.e.,  the  reproduction  process);  and  4)  the 
crossover  and  mutation  processes  are  computed  to  complete  one 
generation  of  GA.  This  process  is  repeated  to  produce  the  species 
with  the  highest  fitness  value,  as  shown  in  the  evolutionary 
process  in  Fig.  17.  GA  has  been  used  to  optimize  the  parameters 
of  other  algorithms  (such  as  FLC  or  ANNs);  only  Shaiek  et  al.  [49] 
used  GA  to  solve  PSC  problems.  Compared  with  other  EC  methods, 
implementing  GA  is  complex  and  difficult  to  achieve  using  low- 
cost  microcontrollers.  Consequently,  Ref.  [49]  only  provided 
simulated  results.  Ref.  [50]  is  one  conference  paper  that  used  GA 
concepts. 

3A.3.4.  Differential  evolution.  DE  is  also  an  EC  method,  and  uses 
architecture  similar  to  that  of  GA.  The  primary  difference  is  that  DE 
uses  special  differential  operators  to  replace  the  crossover  process 
in  GA  and  produce  the  next  generation.  Comparing  to  GA,  DE  is 
simple  to  implement  and  requiring  few  parameter  adjustments. 
Currently,  no  journal  paper  has  used  DE  to  solve  PSC  problems. 


Fig.  16.  Basic  FLC  framework. 


Fig.  17.  The  evolutionary  process  of  GA. 
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Ref.  51]  and  Ref.  [52]  are  conference  papers  based  on  DE,  and  they 
provided  only  simulated  results. 


3.1. 3. 5.  Ant  colony  system.  An  ACS  is  an  EC  method  based  on  swarm 
intelligence;  its  primary  advantage  is  immediately  adapting  command 
values  according  to  environmental  changes.  Thus,  ACSs  are  suitable  for 
conducting  MPP  tracking  in  changing  environments.  In  an  ACS,  each 
agent  selects  its  path  randomly  at  first.  If  the  path  the  agent  chooses  is 
short  (has  high  fitness  value),  the  agent  leaves  concentrated 
pheromone  on  the  path.  In  the  next  iteration,  the  agent  chooses  its 
path  based  on  the  concentration  of  pheromone  on  that  path.  The  more 
concentrated  the  pheromone  is,  the  higher  probability  is  that  the 
agent  chooses  that  path.  The  probability  of  the  agent  moving  from 
position  i  to  position  j  is  as  follows: 


2 

s  e  allowedk 


(13) 


where  t,-j  is  the  pheromone  quantity  on  path  (ij),  rpj  is  the 
initial  value  of  pheromone  on  path  (ij),  a  and  (3  are  adjustable 
parameters,  and  the  updated  pheromone  equation  can  be 
expressed  as  follows: 

Tij(t+n)=piij(t)+ATij  ,  Atij=  Y,  Art  (14) 

k  =  1 


In  (14),  At\-  represents  the  quantity  of  pheromone  that  the  kth 
agent  left  on  path  (ij)  in  this  iteration,  whereas  rjj  represents  the 
total  pheromone  quantity  left  on  path  (ij).  p  is  a  coefficient  such 
that  (1  -  p)  represents  the  evaporation  rate.  The  frequently  used 
At\-  calculation  method  is  as  follows: 

k  (  &  ,  if  the  k  -thant  passes  path  ( ij )  in  this  iteration  ) 

=  1  0  ,  other  | 


(15) 


In  (15),  Q.  is  a  constant  and  Lk  represents  the  distance  of  the 
path  that  the  kth  agent  passes  in  this  iteration. 

Based  on  the  principle  represented  by  the  mentioned  equa¬ 
tions,  Jianga  et  al.  [53]  used  ACS  to  solve  GMPPT  problems. 
However,  ACSs  are  difficult  to  implement  using  low-cost  micro¬ 
controllers;  therefore,  Ref.  [53  provided  only  simulated  results. 

3. 1.3. 6.  Particle  swarm  optimization.  PSO  is  also  an  EC  method 
based  on  swarm  intelligence.  In  the  PSO  method,  each  particle  is 
defined  by  its  own  position  and  velocity.  The  behavior  of  particles 
within  the  swarm  is  influenced  by  the  experiences  of  neighboring 
particles.  Each  particle  follows  the  current  best-performing  particle  to 
search  within  the  solution  space.  In  the  PSO  method,  particles  are 
initially  randomly  or  evenly  distributed  in  the  solution  space.  In  each 
iteration,  particles  follow  two  “extreme  values”  to  update  themselves. 
The  first  value  is  the  optimal  solution  that  the  particle  found,  also 
called  the  personal  best  value  (pbest);  the  second  extreme  value  is  the 
optimal  solution  discovered  by  the  entire  swarm,  also  called  the  global 
best  value  (gbest).  After  these  optimal  solutions  are  found,  the 
particles  update  their  velocity  and  position  based  on  the  following 
equations,  as  shown  in  Fig.  18: 

v|<+ 1  =wv]f  +  c-lrand\(pbestl- -x'-l  +  Cjrand^igbestb -xf)  (16) 


7 


learning  factors),  which  were  used  to  direct  the  search  toward  the 
local  and  global  optimal  solutions,  respectively;  and  randA  and 
rand2  are  random  numbers  between  [0, 1  ]. 

Of  all  the  EC  methods,  the  PSO  algorithm  is  most  widely 
applied  to  GMPPT.  Because  PSO  is  an  optimization  method  based 
on  swarming,  it  can  conduct  GMPPT  in  distributed  architecture 
[54-56]  or  centralized  architecture  [57-60].  Regarding  algorithms, 
Chowdhury  and  Saha  [54]  used  adaptive  perceptive  PSO  (APPSO), 
whereas  the  remaining  references  used  basic  PSO  algorithms. 
Because  implementing  APPSO  is  complex,  Ref.  [54]  provided  only 
simulated  results,  whereas  the  remaining  references  provided 
experimental  results.  This  observation  shows  that  the  PSO  algo¬ 
rithm  is  an  EC  method  that  can  be  readily  implemented  using 
microcontrollers.  In  terms  of  centralized  architecture,  the  PSO  method 
can  accurately  track  the  GMPP  in  various  shading  patterns,  conduct 
direct  cycle  control  [57-60 ],  and  attain  reduced  steady-state  oscillation 
[58].  In  [60],  random  number  in  the  accelerations  coefficient  is 
removed,  developing  a  deterministic  PSO  (DPSO)  mechanism  that 
improved  the  tracking  speed.  In  addition,  Liu  et  al.  [59]  proposed 
suggestions  regarding  practically  implementing  PSO. 

Conference  papers  [61-64]  involved  similar  PSO  methods; 
authors  in  [63]  and  [64]  combined  PSO  and  Bayesian  fusion, 
rendering  a  system  that  can  yield  excellent  GMPPT  results  in  a 
changing  environment. 

3. 1.3. 7.  Chaos  search  method.  CS  is  a  stochastic  search  method 
based  on  chaos  theory,  using  the  ergodicity  of  the  chaos  variable 
to  conduct  an  optimization  search.  Consequently,  the  search 
results  of  this  method  are  superior  to  those  of  other  search 
methods  using  pure  random  numbers.  Zhou  et  al.  [65]  used  CS 
method  based  on  dual  carrier  to  conduct  GMPP  search,  employing 
(18)  and  (19)  to  randomly  produce  control  variables,  measure 
output  power,  and  find  the  GMPP: 

x- =  a-\-xn(b-a)  (18) 

yri  =  a  +j(yn  +  2  )(b-a)  (19) 

According  to  [65],  CS  can  accurately  find  the  GMPP  under  PSC; 
however,  simulated  and  experimental  results  did  not  verify  this 
theory. 


xk+\  =xk  +  vk  +  \ 


(17) 


Where  v\  is  the  velocity  of  particle  i  of  the  kth  iteration;  w  is  the 
weight  value;  cA  and  c2  are  acceleration  coefficients  (also  called 


3.I.3.8.  Tabu  search  method.  The  Tabu  search  method  used  records 
of  previous  search  results  to  avoid  falling  into  local  optimal 
solutions.  Zheng  et  al.  [66]  proposed  using  tabu  search  to  solve 
GMPP  search  problems;  however,  the  authors  did  not  conduct  an 
in-depth  exploration  of  its  implementation  or  design  guidelines. 
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3.1.4.  Other  methods 

This  subsection  introduces  additional  GMPPT  techniques  that 
were  not  categorized  into  the  three  previous  subsections. 

3. 1.4.1.  Extreme  seeking  control  method.  The  extreme  seeking 
control  (ESC)  method  is  an  adaptive  close-loop  control  method 
used  to  search  for  extreme  values.  Compared  with  traditional  P&O, 
ESC  enhances  tracking  speed.  Fig.  19  shows  the  block  diagram  of  a 
typical  ESC.  To  address  PSC,  Lei  et  al.  [67]  used  a  method  similar  to  that 
of  [30]  to  sweep  and  search  for  the  GMPP  location.  The  largest 
difference  between  this  method  and  that  of  [30]  is  that  it  replaced 
P&O  by  using  ESC.  Heydari-doostabad  et  al.  68]  improved  the  tracking 
speed  of  ESC  and  the  startup  power  drop,  employing  a  method  similar 
to  that  of  [28]  to  rapidly  move  the  OP  close  to  the  GMPP.  This  method 
exhibits  satisfactory  tracking  speed;  however,  ESC  calculation  is 
complex,  the  internal  high-pass  filter  is  susceptible  to  noise 
interference,  and  actual  implementation  is  difficult.  In  addition,  ESC 
replaced  P&O  to  provide  the  control  signal  required  for  the  power 
switch,  hindering  ESC  integration  into  traditional  PGS  firmware. 
Consequently,  this  method  is  limited  to  simulation  in  current 
literature,  and  actual  experimental  data  cannot  be  provided. 

3. 1.4.2.  Model  predictive  control  method.  The  model  predictive 
control  (MPC)  method  can  be  used  to  forecast  future  events  and 
implement  controls  based  on  those  forecasts;  therefore,  MPC 
involves  optimizing  based  on  the  control  quantity  required  in 
the  current  situation.  Bouilouta  et  al.  [69]  combined  MPC  and 
large-step  GMPP  search  to  address  PSC  problems.  The  MPC 


method  exhibits  favorable  tracking  speed  and  accuracy,  although 
complex  calculations  are  required.  Furthermore,  the  MPC  design 
parameters  are  related  to  the  system  model  and  the  power 
converter  topology.  Compared  with  MPC,  the  design  of  other 
methods  need  not  take  the  power  converter  topology  into 
account;  thus,  actual  implementation  of  MPC  is  system  dependent 
and  difficult.  Consequently,  work  in  [69]  only  provided  simulated 
results. 


3.1.5.  Comparing  firmware-based  methods 

In  this  subsection,  a  simple  comparison  will  be  conducted, 
incorporating  the  journal  papers  that  used  the  aforementioned 
firmware-based  methods  suitable  for  use  in  centralized  architec¬ 
ture.  The  comparison  addresses  the  dependence  on  system  archi¬ 
tecture,  tracking  speed,  probability  of  successfully  finding  the 
GMPP,  complexity  of  implementation,  and  sensors  required. 
Table  1  summarizes  the  comparison  results,  providing  a  reference 
for  PGS-designers  and  researchers  during  method  selection. 


3.1.6.  Methods  for  detecting  partial  shading  condition 

The  method  for  detecting  PSC  is  often  overlooked  in  GMPPT- 
related  literature  mainly  because  the  transient  state  of  various 
operating  conditions  need  not  be  considered  (e.g.,  operating 
conditions  changing  from  uniform  insolation  to  PSC  or  vice  versa) 
when  writing  a  paper.  However,  PGS  design  engineers  must  be 
able  to  detect  whether  the  current  system  is  exposed  to  uniform 
insolation  or  PSC,  and  whether  the  system  operating  condition 
changed.  Because  MPPT  techniques  for  handling  PSC  problems 
must  address  numerous  factors,  such  techniques  are  typically 
slower  compared  with  traditional,  rapid  MPPT  methods  (e.g.,  the 
variable-step  P&O  method).  Therefore,  an  accurate  understanding 
of  the  current  operating  status  facilitates  selecting  suitable  algo¬ 
rithms  and  avoiding  extending  the  tracking  time.  In  this  subsec¬ 
tion,  the  PSC  detection  methods  proposed  by  the  references  listed 
in  Table  1  are  compared.  Table  2  summarizes  the  comparison 
results,  showing  whether  the  proposed  detection  methods  can 
accurately  distinguish  between  uniform  insolation  and  PSC,  and 
whether  they  can  determine  if  the  GMPP  changed  and  subse¬ 
quently  reinitialize  tracking  for  the  five  possible  operating  condi¬ 
tion  changes.  Fig.  20  shows  the  five  possible  changes  in  operating 
condition  tested  in  this  study. 


Table  1 

Comparison  of  firmware-based  methods. 


Method 

Dependence 
on  system 

Tracking  speed 

Hit  count 

Algorithm  complexity 

Sensors  required 

Notes 

[26-28] 

Y 

Fast 

Low 

Low 

V,I 

[30] 

N 

Slow 

Medium 

Medium 

V,I 

[31] 

N 

Medium 

Medium  to  high 

Low 

V,I 

Can  easily  be  integrated  to  original  PGS 

[32], [33] 

N 

(relate  to  alpha) 
Fast 

Medium  to  high 

Low 

V,I 

Can  easily  be  integrated  to  original  PGS 

[36] 

N 

Fast 

High 

Medium 

V,I 

Requires  recursive  operation 

[37], [38] 

N 

Fast 

High 

Medium 

V,I 

Requires  recursive  operation 

[43], 

Y 

Fast 

Relate  to  the  completeness 

Low  to  medium 

Irradiance 

[44] 

[47] 

N 

Fast 

of  the  training  data 

Medium  to  high 

Medium 

Temperature 

V,I 

[49] 

N 

Slow  to  medium 

High 

Medium  to  high 

V,I 

Stochastic  process 

[53] 

N 

Slow  to  medium 

High 

Medium  to  high 

V,I 

Stochastic  process 

[57-60] 

N 

Slow  to  medium 

High 

Medium 

V,I 

Stochastic  process 

[67] 

N 

Slow  to  medium 

High 

Medium  to  high 

V,I 

Stochastic  process 

Table  2 

Comparing  PSC  detection  methods. 


Method 


Ref.  Parameter  description 


Case  I 


Case  II 


Case  III 


Case  IV 


Case  V 


Reinitial  Uniform/Non-  Reinitial  Uniform/Non-  Reinitial  Uniform/Non-  Reinitial  Uniform/Non-  Reinitial  Uniform/Non¬ 
uniform  uniform  uniform  uniform  uniform 


/(Sfc+1)-J(S,<)^ 

KSk) 


P(/<+l)-P(/<) 

P(k) 


V(k+l)-V(k)<AV 

P(fc+1)-P(fe) 

P(/<) 


V(k+ 1)|  <  -AV 

P(fe+1)-P(fe) 

P(k) 


[53] 

J(S k+i):  The  current  of  the  (/<+l)-th 
iteration 

Yes 

No 

Yes 

No 

Yes 

No 

Yes 

No 

No 

No 

[59] 

/(5k):  The  current  of  the  k-th  iteration 

A/:  A  pre-defined  constant 

P(/c+l):  The  power  of  the  (/<+l)-th 
iteration 

Yes 

No 

Yes 

No 

Yes 

No 

Yes 

No 

No 

No 

[49] 

P(/<):  The  power  of  the  k-th  iteration 

AP:  A  pre-defined  constant 

V(/c+l),  P(k+1):  The  voltage/power  of 
the  (/<+l)-th  iteration 

Yes 

No 

Yes 

No 

Yes 

No 

Yes 

No 

No 

No 

[56] 

V(k),  P(k ):  The  voltage/power  of  the  k-th 
iteration 

AV,  AP:  Pre-defined  constants 

V(k+ 1),  P(k+1):  The  voltage/power  of 

Yes 

No 

Yes 

No 

Yes 

No 

Yes 

No 

No 

No 

the  (fc+l)-th  iteration 


>0.1 

ld3 


Vd!  -  Vd2 

Vdi 


>0.2 


P'l  (X2  d )  =  ^X2~d^~^>k~'l^x'l+d^ 


Pfc(X2-d)-Pfc-l(Xl+d) 

P'k(X2-d) 


P(k):  The  power  of  the  k-th  iteration 
AP:  Pre-defined  constants 

Vdl,Jdl:  The  voltage/current  sample  of  Yes  Yes  Yes  Yes  Yes  Yes  Yes  Yes  No  No 

point  1  (in  the  middle) 


[57]  Vd2,Id2:  The  voltage/current  sample  of 
point  2  (at  the  left  side) 

[60]  Vd3,Id3:  The  voltage/current  sample  of 
point  3  (at  the  right  side) 

[31,37]  d:  The  searching  direction  Yes  Yes  Yes  Yes  Yes  Yes  Yes  Yes  No  No 


Pk(*2_d):  New  obtained  power 
Pk_i(x1  +  d):  Previous  obtained  power 
P'k(*2-d):  Averaged  power 
r:  A  pre-defined  constant 

AP>APcrit  [30,33]  AP:  The  variation  of  the  power  Yes  No  Yes  No  Yes  No  Yes  No  No  No 
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3.2.  Hardware-based  GMPPT  techniques 

In  addition  to  using  firmware-based  methods  to  address  the 
multiple  peak  values  caused  by  PSCs,  many  researchers  have 
proposed  using  hardware  methods  to  solve  this  problem.  These 
methods  are  simply  discussed  as  follows. 


3.2.2.  Investigation  on  system  configuration  and  reconfiguration 
technique 

If  unmanaged,  PSC  cause  the  useable  power  of  a  system  to 
decrease.  Works  in  [70-72]  studied  the  three  configuration  meth¬ 
ods  used  in  large  PGSs,  including  series-parallel  (SP),  bridge-linked 
(BL),  and  total-crossed-tied  (TCT).  The  results  showed  that  in  PSC, 
TCT  and  BL  attain  superior  performance  compared  with  SP. 
Comparing  TCT  and  BL,  TCT  produces  more  power.  Furthermore, 
works  in  [73]  and  [74]  proposed  immediately  adjusting  the 
configuration  of  PV  arrays  according  to  the  shading  pattern  to 
compensate  for  the  power  losses  caused  by  PSC.  Although  recon¬ 
figuration  method  can  mitigate  the  effects  of  PSC,  this  technique 
requires  using  a  switch  matrix  to  implement  architecture  changes; 
therefore,  the  system  is  expensive  and  the  controller  design  is 
complex.  In  addition,  reconfiguration  technique  cannot  guarantee 
the  successful  tracking  of  GMPP  in  all  shading  patterns.  Ref.  [75] 
and  Ref.  76]  are  conference  papers  discussed  system  configura¬ 
tions,  whereas  Ref.  [77-81]  are  conference  papers  explored  using 
the  reconfiguration  method  to  reduce  the  effects  of  PSC. 
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3.2.2.  Multi-level  converter 

When  a  system  adopts  centralized  architecture  as  shown  in 
Fig.  4(a),  a  multi-level  converter  can  be  used  when  each  PV 
module  in  the  PGS  must  be  independently  controlled  [82-84]. 
Multi-level  converter  can  address  PSC  problems;  however,  because 
its  number  of  levels  is  proportional  to  the  number  of  series- 
connected  PV  modules  in  a  PGS,  after  the  architecture  of  PGS 
changes,  the  entire  hardware  system  and  controller  must  be  re¬ 
designed.  In  addition,  the  greater  the  number  of  series-connected 
modules  is,  the  more  power  semiconductor  devices  are  required, 
increasing  costs  and  power  losses.  Ref.  [85-89]  are  conference 
papers  explored  multi-level  converters. 


3.2.3.  Distributed  architecture 

Fig.  4(b)-(d)  show  that  in  addition  to  using  a  centralized 
architecture,  a  PGS  can  be  designed  to  apply  a  DC/DC  converter 
or  DC/AC  inverter  for  each  PV  module.  This  architecture  is  called  a 
distributed  architecture.  Because  each  power  circuit  has  its  own 
designated  controller,  even  if  PV  modules  in  the  system  experience 
various  levels  of  shading,  the  distributed  architecture  can  continue 
to  track  the  GMPP  of  the  entire  system.  Typically,  when  a  DC/DC 
converter  is  attached  to  the  back  of  the  PV  module  as  shown  in 
Fig.  4(b)  and  (c),  this  architecture  is  called  a  module-integrated 
converter  (MIC)  [90-96].  Conversely,  when  a  DC/AC  inverter  is 
attached  to  the  back  of  the  PV  module  as  shown  in  Fig.  4(d),  this 
architecture  is  often  called  a  micro-inverter  [97  .  Distributed 
architecture  can  guarantee  tracking  of  the  system  GMPP;  further¬ 
more,  because  its  control  is  distributed,  the  reliability  and  robust¬ 
ness  of  the  system  increase.  Flowever,  the  hardware  costs  of 
distributed  architecture  are  often  high,  and  implementing  the 
controller  is  complex  and  requires  a  high  processing  speed 
microcontroller.  In  addition,  because  micro-inverters  are  directly 
connected  to  the  utility,  the  hardware  design  must  address 
grid-interconnection  and  safety  problems  [97].  Ref.  [98-104]  are 
conference  papers  studied  MICs,  whereas  Ref.  105-107]  are 
conference  papers  explored  micro-inverters. 
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Case  I 


Case  II 


Case  III 


Case  IV 


Case  V 


I  I 


\  \  I 


Fig.  20.  The  five  possible  changes  in  operating  mode  tested  in  this  study. 


hub  ^mod 


Fig.  21.  Equalizer-assisted  architecture. 

3.2.4.  Equalizer- assisted  architecture 

Equalizer-assisted  architecture  as  shown  in  Fig.  21  has  been  a 
popular  method  of  addressing  PSC  problems  in  the  past  two  years. 
Fig.  21  shows  that  the  operating  principle  of  this  method  is  similar 
to  the  parallel-connected  MIC  shown  in  Fig.  4(b).  The  primary 
difference  is  that  parallel-connected  MIC  must  process  all  the 
power  supplied  by  the  PV  module,  whereas  the  equalizer-assisted 
architecture  must  only  process  the  power  difference  between 
unbalanced  PV  modules.  Therefore,  this  architecture  guarantees 
tracking  of  the  system  GMPP  and  its  hardware  costs  are  less 
compared  with  that  of  the  parallel-connected  MIC  method  [108- 
112  .  Because  this  method  is  similar  to  the  equalizer  used  in  the 
series-connected  batteries  of  electric  vehicles,  researchers  have 
called  this  architecture  a  PGS  equalizer.  The  advantages  and 
disadvantages  of  this  method  are  the  same  as  those  of  distributed 
architecture.  Conference  paper  [113]  also  examined  equalizer- 
assisted  architecture. 


3.2.5.  Summary  and  comparison  of  hardware  methods 

Previous  descriptions  show  that  numerous  types  of  hard¬ 
ware  architecture  can  solve  or  mitigate  PSC  problems.  Because 
hardware-based  solutions  relate  to  the  selection  of  power  con¬ 
verters,  which  must  be  designed  to  address  system  specifications 
and  power  ratings,  and  most  system  specifications  used  in 
literature  are  distinct.  Therefore,  no  further  in-depth  discussions 
are  conducted  regarding  hardware-based  solutions;  however, 
relevant  information  regarding  such  solutions  is  summarized  in 
Table  3. 


4.  Simulated  results  of  selected  methods 

The  number  of  references  and  the  number  of  citations  demon¬ 
strated  that  the  most  commonly  used  firmware-based  techniques 
for  addressing  PSC  were  the  two-stage,  segmental  search,  and  SC 
methods.  In  this  subsection,  five  representative  methods  were 
selected  and  MATLAB  was  employed  to  test  these  methods  and 
perform  comparisons.  The  methods  selected  in  this  study  were 
Ref.  [28]  (representing  the  system  characteristic  curve  method, 
cited  73  times),  Ref.  30]  (the  most  cited  method,  cited  215  times), 
Ref.  [32]  (representing  the  two-stage  searching  method,  cited  19 
times),  Ref.  [36]  (representing  the  segmental  search  method,  cited 
55  times),  and  Ref.  [58]  (representing  the  SC  method,  cited  45 
times).  Because  the  parameter  settings  of  30]  and  [32]  influenced 
the  GMPPT  performance,  simulations  of  [30]  and  [32]  were 
conducted  using  different  parameter  settings  to  assess  the  tracking 
performance.  Table  4  lists  the  various  methods  and  relevant 
parameter  setting  values.  In  Section  2,  it  was  demonstrated  that 
distinct  shading  patterns  produce  distinct  P-V  characteristic 
curves.  However,  most  references  tested  only  a  few  specific 
shading  patterns;  thus,  the  proposed  methods  cannot  be  guaran¬ 
teed  suitable  for  all  possible  shading  patterns.  In  this  study, 
simulations  were  performed  based  on  a  5slp  system  constituting 
five  Sun  Power  E19/240  PV  modules,  and  MATLAB  was  used  to 
simulate  100-1000  W/m2  of  irradiance  levels  on  each  module  in 
100  W/m2  increments.  Consequently,  each  module  had  10  possible 
irradiance  levels.  To  produce  a  P-V  characteristic  curve  that 
exhibited  five  peak  values  in  a  5slp  system,  the  irradiance  level 
for  each  module  must  be  distinct.  In  five  series-connected  modules 
that  exhibit  different  irradiance  levels,  there  are  C(  10,5)  =  252 
possible  shading  patterns.  Therefore,  the  eight  methods  listed  in 
Table  4  were  tested  using  all  252  possible  shading  patterns  in  this 
paper.  Table  5  lists  the  results.  The  specification  of  the  PV  module 
used  in  this  study  is  summarized  in  Table  6. 
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Table  3 

Comparison  of  hardware-based  methods. 


Method 

Advantage 

Disadvantage 

Investigation  on  system  configuration  [70-72] 
Reconfiguration  methods  [73-74] 

No  hardware  required 

No  power  converter  required 

Can  only  mitigate  PSC  effect 

Cannot  guarantee  GMPP  tracking 

High  system  cost 

Complicated  control 

Multi-level  converter  [82-84] 

MPPT  for  individual  PV  module 

Individual  PV  modules  may  not  operate  at  their  MPP. 
System  should  be  re-designed  when  PGS  architecture 
changes 

High  switching  loss 

Distributed  power  structure  (series/parallel  MIC) 

[90-96] 

MPPT  for  individual  PV  module 

High  GMPPT  accuracy 

System  upgradable 

Higher  hardware  costs 

Complicated  control 

Distributed  power  structure  (micro-inverter)  [97] 

MPPT  for  individual  PV  module 

High  GMPPT  accuracy 

System  upgradable 

Compact  size 

Higher  hardware  costs 

Complicated  control 

Grid-interconnection  issues  should  be  taken  into 

account 

Balancer  auxiliary  structure  [108-112] 

MPPT  for  individual  PV  module 

High  GMPPT  accuracy 

Only  process  the  power  difference  between  mismatched  PV 
modules 

System  upgradable 

Compact  size 

Higher  hardware  costs 

Complicated  control 

Table  4 

The  methods  simulated  in  this  study  and  their  corresponding  parameter  settings. 


Method 

Ref. 

Initial  Settings 

Method  1 

36 

P&O  perturbation  step  in  2nd  stage:  1.0  V 

Method  2 

30 

Search  step  in  1st  stage:  0.6  Voc,one 

P&O  perturbation  step  in  2nd  stage:  1.0  V 

Method  3 

30 

Search  step  in  1st  stage:  0.7  Vocone 

P&O  perturbation  step  in  2nd  stage:  1.0  V 

Method  4 

30 

Search  step  in  1st  stage:  0.8  Voc>one 

P&O  perturbation  step  in  2nd  stage:  1.0  V 

Method  5 

32 

Search  step  in  1st  stage:  20  W 

P&O  perturbation  step  in  2nd  stage:  1.0  V 

Method  6 

32 

Search  step  in  1st  stage:  50  W 

P&O  perturbation  step  in  2nd  stage:  1.0  V 

Method  7 

28 

N.A. 

Method  8 

59 

ci)=0.4,  Ci  =  1.2,  c2  =  1.6 

In  Table  5,  hit_count  refers  to  the  number  of  successful  times  the 
GMPP  was  located  during  the  252  tests.  Because  the  second  stages  of 
[28],  [30],  and  [32]  used  P&O  (locating  MPP  using  P&O  yields 
oscillations),  tracking  success  was  defined  as  the  located  MPP  being 
99.8%  larger  than  the  actual  GMPP;  otherwise,  the  result  was  regarded 
as  tracking  failure.  In  certain  shading  patterns,  the  difference  between 
the  highest  and  second  highest  peaks  in  the  P-V  characteristic  curve 
was  small  (such  as  shading  pattern  [800  W/m2, 400  W/m2, 300  W/m2, 
200  W/m2, 100  W/m2],  GMPP= 203.9  W,  second  highest  peak=  195.0 
W);  this  can  cause  errors  in  certain  GMPPT  methods.  However,  the 
difference  in  the  obtained  power  will  be  small  if  the  system  tracks  the 
second  highest  peaks  instead  of  the  highest  ones.  Consequently, 
obtaining  results  that  presented  a  low  hit_count  and  a  relatively  high 
output  power  was  likely.  Therefore,  tracking  accuracy  should  also  be 
taken  into  account  for  performance  evaluation.  Table  5  lists  a 
comparison  of  tracking  accuracy  for  the  eight  methods.  In  Table  5, 
avg_acc  refers  to  the  average  accuracy  value  of  MPP  tracking  for  the 
252  tests;  that  is,  Jn  =  -iMPPacc,n/ 252,  and  MPPaccn  is  defined  as  the 
ratio  of  the  MPP  found  by  the  algorithm  divided  by  the  actual  GMPP  in 
the  nth  test.  Similarly,  max_acc  and  min_acc  refer  to  the  maximal  and 


minimal  values  of  the  252  test  results,  respectively.  In  addition  to 
tracking  accuracy,  tracking  speed  is  also  critical  to  system  design. 
When  comparing  tracking  speeds  of  existing  methods,  many  research¬ 
ers  have  compared  actual  tracking  times  (unit  in  seconds);  however, 
unless  the  PGS  specifications  and  power  converter  of  two  systems  are 
identical  (indicating  that  the  hardware  settling  time  is  consistent),  this 
comparison  provides  no  value.  For  example,  the  practical  design  of  an 
MPPT  system  must  address  the  settling  time  of  power  circuits  to  select 
the  interval  between  two  sampling  times.  If  the  sampling  speed  is 
faster  compared  with  the  system  settling  time,  the  value  of  the  voltage 
and  current  sampled  by  the  algorithm  does  not  indicate  steady  state, 
and  thus  cannot  be  used  to  assess  the  effects  generated  by  the  control 
command.  However,  the  settling  time  of  a  low-power  switching 
converter  is  shorter  than  that  of  a  high-power  one.  Therefore,  despite 
using  identical  GMPPT  algorithms,  the  tracking  time  of  a  low-power 
system  may  be  shorter  compared  with  that  of  a  high-power  system.  In 
other  words,  assessing  the  tracking  speed  of  MPPT  methods  by 
comparing  tracking  times  but  without  considering  other  factors 
produces  biased  results.  Therefore,  the  authors  of  this  paper  propose 
that  the  number  of  iterations  required  (related  to  only  the  method) 
rather  than  the  time  (related  to  the  system  and  the  method)  should  be 
used  to  compare  tracking  speeds.  Iteration  refers  to  the  number  of 
times  a  specific  function  is  calculated.  Using  the  P&O  method  as  an 
example,  one  iteration  refers  to  one  perturbation.  In  another  example, 
the  number  of  iterations  is  the  number  of  MPP  calculation  interrupts 
when  a  proposed  method  is  implemented  by  microcontrollers.  In 
Table  5,  avg_step  refers  to  the  average  number  of  required  iterations 
for  tracking  the  MPP  in  the  252  tests.  Similarly,  max_step  and 
min_step  refer  to  the  maximal  and  minimal  required  iteration 
numbers  in  the  252  MPPT  tests,  respectively.  Of  the  eight  testing 
methods,  only  Method  8  (PSO  algorithm)  exhibited  stochastic  char¬ 
acteristics;  that  is,  the  tracking  results  obtained  in  each  execution 
could  vary.  Thus,  test  results  of  Method  8  provided  in  Table  5  is  the 
average  of  20  repeated  tests  using  the  same  shading  pattern;  conse¬ 
quently,  the  max_step  and  min_step  of  Method  8  are  not  integer 
numbers. 

From  Table  5,  the  following  can  be  concluded: 

1)  P&O,  IncCON,  and  variable-step  P&O/IncCON  are  the  most 
commonly  used  methods  in  the  second  stage  of  the  two- 
stage  method.  Therefore,  the  performance  of  the  two-stage 
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Table  5 

Simulated  results. 


Direct 

Patel_0.6 

Patel_0.7 

Patel_0.8 

Power_20 

Power_50 

Real 

PSO 

avg_MPPT_acc  (%) 

98.99 

51.36 

51.35 

98.41 

99.85 

99.45 

52.09 

99.96 

max_MPPT_acc  (%) 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

99.79 

100.00 

min_MPPT_acc  (%) 

92.48 

17.43 

17.43 

68.51 

96.69 

81.21 

31.45 

99.20 

hit_count 

183 

11 

11 

153 

227 

215 

1 

251.85 

avg_MPPT_step 

20.94 

71.44 

71.87 

473.15 

33.98 

23.58 

42.26 

129.19 

max_MPPT_step 

31.00 

102.00 

112.00 

1185.00 

53.00 

38.00 

47.00 

700.00 

min_MPPT_step 

16.00 

60.00 

65.00 

38.00 

19.00 

13.00 

22.00 

7.00 

Table  6 

The  specification  of  the  sun  power  E19/240  PV 
module. 


Parameters  for  model 

Maximum  power  (Pm ax) 

240  W 

Open  circuit  voltage  (Voc) 

48.6  V 

MPP  voltage  (Vmp) 

40.5  V 

Short  circuit  current  (/sc) 

6.30  A 

MPP  current  (Jmp) 

5.93  A 

Temperature  coefficient  (av) 

- 132.5  mV/K 

method  is  determined  mainly  based  on  the  first-stage  (GMPP 
location  search). 

2)  The  probability  of  finding  the  GMPP  by  using  a  system 
characteristic  curve  method,  which  is  based  on  PV  panel 
equivalent  resistance,  is  relatively  low.  In  [28],  when  PSC 
yielded  a  complex  P-V  characteristic  curve  (increasing  the 
number  of  LMPPs),  the  probability  of  this  method  finding  the 
GMPP  significantly  decreased  (1/252).  Thus,  this  method 
should  only  be  used  for  simple  P-V  characteristic  curves  (when 
used  to  assess  a  2slp  system,  this  method  can  yield  a  20/45 
probability  of  success  with  avg_acc=85.42%). 

3)  The  method  used  in  [30]  attained  the  highest  number  of  citations 
in  current  GMPPT  paper.  However,  this  method  is  designed  after 
observing  the  P-V  characteristic  curve  caused  by  PSC  and  convert¬ 
ing  such  observations  into  simplified  search  rules.  Thus,  the 
probability  of  finding  the  GMPP  is  also  low.  The  performance  and 
design  considerations  of  this  method  are  outlined  as  follows:  (1)  in 
[30],  the  author  recommended  using  a  0.6-0.7  Voc,one  step  to 
conduct  jump  search.  However,  extensively  testing  by  using  a  5slp 
system  assembled  from  Sun  Power  El 9/240  PV  modules  showed 
that  0.8  Voc,one  yielded  better  results.  Therefore,  prior  to  using  this 
method,  simple  analysis  on  system  characteristics  or  parameter 
setting  experiments  should  be  conducted  to  determine  the  optimal 
design  parameters;  (2)  Although  the  hit_count  (153/252)  for  this 
method  was  not  high,  the  avg_acc  result  (98.41%)  indicated  a  high 
level  of  tracking  accuracy,  suggesting  that  this  method  did  not  track 
the  real  GMPP,  but  the  power  of  the  tracked  LMPP  was  not 
substantially  different  compared  with  the  power  of  the  GMPP.  In 
other  words,  from  the  perspective  of  electricity  sales  profit,  this 
method  is  worth  adopting;  and  (3)  Table  5  shows  that  large 
numbers  of  steps  are  required  in  this  method  (473.15),  indicating 
that  the  tracking  time  is  relatively  long.  Because  this  method  uses 
P&O  to  track  each  LMPP,  improving  this  procedure  can  enhance 
tracking  speed. 

4)  The  method  in  [32]  is  another  representative  method  of  the 
two-stage  searching  method.  The  simulated  results  showed 
that  this  method  attained  excellent  hit  counts,  tracking  accu¬ 
racy,  and  tracking  time.  The  primary  difference  between  this 
method  and  that  used  in  [30]  was  that  this  method  involves 
only  a  brief  scan  of  the  P-V  characteristic  curve  in  the  first 
stage;  therefore,  the  slow  tracking  speed  encountered  in  [30] 
(using  P&O  to  track  each  LMPP  reduces  tracking  speed)  can  be 


avoided.  The  simulated  results  of  this  method  indicated  that 
when  large  steps  were  selected  for  the  sweep  (AP=50  W),  the 
system  required  decreased  time  to  find  the  GMPP  (23.58  steps), 
but  could  miss  the  GMPP  (hit_count=215).  Conversely,  when 
small  steps  were  selected  (AP=20W),  the  system  had  a  high 
probability  of  finding  the  GMPP  (hit_count=227),  but  required 
a  relatively  long  time  (33.98  steps).  Therefore,  using  this 
method  requires  addressing  the  tradeoff  between  step  size 
and  search  time.  Because  this  method  is  simple  to  implement 
and  can  be  readily  integrated  into  original  PGS  firmware,  it  is  a 
method  worth  using,  implying  two  directions  of  future  study. 
First,  finding  a  sweeping  rule  that  can  improve  hit_count,  and 
second,  improving  the  P&O  method  used  in  the  second  stage  to 
enhance  the  GMPPT  speed  (such  as  using  variable-step  P&O). 

5)  The  method  in  [36]  can  also  be  used  to  balance  hit  count,  tracking 
accuracy,  and  tracking  time;  however,  this  method  cannot  be  easily 
integrated  with  original  PGS  firmware  (typically  requiring  com¬ 
plete  algorithm  re-write),  making  it  only  suitable  for  use  in  newly 
developed  PGSs.  Although  the  mathematical  principle  of  this 
method  is  relatively  simple  and  comprehensive,  actual  coding 
yields  challenges.  First,  the  jump  range  for  sampling  points  used 
in  this  method  (VAi  -►  Va2-»  VA3,  see  Fig.  13)  is  the  largest  of  the  five 
methods  tested,  indicating  that  the  firmware  designer  must 
consider  the  hardware  settling  time  to  avoid  excessive  stress  on 
the  power  switch.  In  addition,  this  method  involves  using  the 
recursive  method  to  gradually  reduce  the  search  range;  this 
challenges  some  firmware  engineers  who  are  only  accustomed  to 
sequential  programming. 

6)  The  method  in  58]  represents  the  SC  method.  According  to  the 
simulated  results,  the  tracking  accuracy  of  this  method  was 
100%  (hit_count=252).  Because  PSO  can  be  applied  to  pro¬ 
blems  more  complex  than  those  examined  in  this  study  (e.g., 
the  traveling  salesman  problem),  this  result  was  unsurprising. 
Additionally,  the  avg_acc  of  this  method  was  the  highest  of  the 
tested  methods,  indicating  that  this  method  yielded  the  highest 
electricity  sales  profit.  However,  based  on  the  simulated  results, 
this  method  required  a  relatively  long  tracking  time;  therefore 
is  unsuitable  for  use  in  areas  that  exhibit  rapid  shading  pattern 
changes.  Because  PSO  is  a  stochastic  optimization  method,  the 
number  of  iterations  required  to  track  the  GMPP  for  each 
execution  may  be  different;  this  distinguishes  SC  from  the  other 
methods.  Furthermore,  the  parameters  in  (16)  and  (17)  substan¬ 
tially  affect  the  tracking  performance  when  this  method  is  used. 
Consequently,  optimizing  these  parameters  is  a  topic  worth 
studying.  According  to  [59]  and  [60],  SC  can  be  implemented  by 
employing  low-cost  microcontrollers  such  as  TMS320  or  dsPIC33 
series  digital  signal  controllers. 


5.  Conclusion 

Extensive  recent  literature  shows  that  the  GMPPT  methods 
suitable  for  use  in  PSC  remain  a  popular  research  topic.  Current 
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literature  comprises  numerous  algorithms  and  hardware  architec¬ 
tures  related  to  addressing  PSC  problems;  therefore,  PGS  designers 
may  be  challenged  to  select  appropriate  firmware  or  hardware 
architecture.  In  this  study  novel  GMPPT  techniques  from  diverse 
references  were  collected,  classified,  compared,  and  summarized; 
the  advantages  and  disadvantages  of  these  techniques  are  dis¬ 
played  in  tables.  Regarding  firmware-based  techniques,  compar¬ 
isons  were  conducted  and  the  indices  that  must  be  considered  in 
the  design  were  summarized  such  as  dependence  on  system 
architecture,  tracking  speed,  the  probability  of  successfully  track¬ 
ing  the  GMPP,  algorithm  complexity,  and  required  sensors.  The 
advantages  and  disadvantages  of  the  hardware-based  techniques 
were  also  discussed.  Furthermore,  PSC  detection  and  determina¬ 
tion  methods,  which  have  seldom  been  discussed  in  current 
literature,  were  reviewed  and  compared.  Finally,  MATLAB  was 
used  to  extensively  simulate  five  commonly  used  GMPPT  algo¬ 
rithms  and  explore  their  tracking  performance;  recommendations 
were  then  proposed  for  using  these  methods.  Based  on  the 
simulated  results,  topics  worth  studying  were  also  suggested. 
Firmware-based  GMPPT  methods  can  continue  to  be  enhanced, 
particularly  in  terms  of  maintaining  hit  count  and  tracking 
accuracy  while  accounting  for  tracking  speed,  as  well  as  addres¬ 
sing  problems  related  to  uniform  insolation  and  PSC  and  detecting 
the  transient  state  in  an  actual  system.  To  conclude,  this  study 
should  serve  as  a  valuable  reference  for  PGS  researchers  and 
application  engineers. 
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